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Abstract 

We consider a finite-horizon market-making problem faced by a dark pool that executes 
incoming buy and sell orders. The arrival flow of such orders is assumed to be random 
and, for each transaction, the dark pool earns a per-share commission no greater than 
the half bid-ask spread. Throughout the entire period, the main concern is inventory 
risk, which increases as the number of held positions becomes critically small or large. 
The dark pool can control its inventory by choosing the size of the commission for each 
transaction, so to encourage, e.g., buy orders instead of sell orders. Furthermore, it can 
submit lit-pool limit orders, of which execution is uncertain, and market orders, which 
are expensive. In either case, the dark pool risks an information leakage, which we 
model via a fixed penalty for trading in the lit pool. We solve a double-obstacle impulse- 
control problem associated with the optimal management of the inventory, and we show 
that the value function is the unique viscosity solution of the associated system of quasi 
variational inequalities. We explore various numerical examples of the proposed model, 
including one that admits a semi-explicit solution. 

Keywords: Market-making, inventory risk, impulse-control problem, quasi variational 
inequality, viscosity solutions. 


1 Introduction 

Market makers are liquidity providers. They set bid and ask quotes and trade with impatient 
investors who seek to immediately buy or sell a certain quantity of a stock by market orders. 
The market maker is willing to hold a non-zero inventory while earning the spread from 

The authors thank J. Walton for useful discussions at the beginning of this work and are grateful to C.A. 
Garcia Trillos for comments which improved this paper. 
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each transaction. Holding a non-zero inventory carries an intrinsic risk associated to the 
unpredictable changes to which the stock price is subject. This risk is further increased 
by a potential information asymmetry due to which the market maker trades in the wrong 
direction. 

In this paper we consider the optimal management of a dark pool inventory and the 
poofs role as a market maker. A dark pool is an alternative trading venue where partic¬ 
ipants do not release their identity and benefit of advantageous prices. While optimal- 
liquidation problems involving dark pools have been treated in recent works by e.g., Kratz 
& Schoneborn (2013), Horst & Naujokat (2014), Crisafi & Macrina (2014), and Graewe 
et al. (2015), as far as we know dark pools have not been studied in a market-making 
context elsewhere. In our situation the dark pool (i) executes incoming buy and sell or¬ 
ders posted by its clients and (ii) may post limit and market orders in a lit pool to control 
the level of its inventory. This view can be justified by the increasing competition between 
dark platforms and the reputational benefit deriving from the speed of execution and price 
advantage. Also, by executing large orders placed by its (institutional) clients, the dark 
pool acquires “reserved” information about the traded asset. Furthermore, the dark pool 
would rather avoid the situation where an investor (client) resorts to the lit pool and thus 
potentially moves the price against the dark pool’s holding. We formulate a double-obstacle 
impulse-control problem and we use the viscosity theory to characterise the solution of the 
associated system of quasi variational inequalities (QVIs). We refer the reader to Crandall 
et al. (1992), Fleming & Soner (2006) and Pham (2009) for a comprehensive treatment of 
viscosity solutions. 

Previous work on market making includes Amihud & Mendelson (1980) who, based 
on Carman (1976), relate the bid-ask prices to the share holding of a risk-neutral agent. 
They find a relationship between the optimal quotes and the distance from the “preferred” 
inventory position. Stoll (1978) considers a two-period model in which a risk-averse agent 
supplies liquidity and maximises his expected utility. Ho & Stoll (1981) use the dynamic 
programming principle (DPP) to obtain the optimal quotes which maximise the terminal 
wealth in a single-dealer market. The recent evolution in financial markets, arisen with 
algorithmic and high-frequency trading, has shifted the optimal market-making problem 
to an order-driven market environment where optimal quotes and trading strategies are 
computed and submitted by electronic machines. For example Avellaneda & Stoikov (2008) 
adapt the market-making framework by Ho & Stoll (1981) to a limit order book (LOB). 
They consider the maximisation of the agent’s expected terminal wealth and consider both, 
the finite and the infinite-time cases. They model the arrival of buy and sell orders by 
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Poisson processes and the dynamics of the mid-price by an arithmetic Brownian motion. 
They find the HJB PDE by means of the DPP and propose an approximation of the optimal 
quotes via asymptotic expansions. This type of problem has been investigated elsewhere, 
too. The works by Cartea & Jaimungal (2012) on risk metrics and by Cartea et al. (2013, 
2014) consider ambiguity and self-exciting processes, respectively. Gueant et al. (2013) 
deal with the inventory risk and reduce a complex optimisation problem to a system of 
ODEs. Guilbaud & Pham (2013) consider a market maker who continuously submits limit 
orders at the best quoted prices and resorts to market orders when the inventory becomes 
too large. They numerically solve a finite-time impulse-control problem and find the optimal 
order sizes and quotes to be posted in the lit pool. 

The paper is organised as follows. After the introduction, in Section we present the 
inventory and cash processes of the dark pool and we describe its trading strategies. In 
Section!^ the optimisation problem is introduced and, by making use of DPI? we derive the 
HJB equation. Section is devoted to a numerical analysis of the dark pool strategies. All 
propositions, necessary to characterise the value function by the unique viscosity solution of 
the system of QVIs, and their proofs can be found in the appendix. 

2 Dark pool as a market maker 

We consider a dark pool that executes incoming buy and sell orders by its clients over a finite 
period of time 0 < t < T < oo and that may resort to the lit pool if its inventory becomes 
critically small or large. As a reward for the service provided, the dark pool chooses a 
per-share commission 5^ for sell orders and 5^ for buy orders, which shall not exceed the 
lit-pool half spread. The dark-pool order flow is affected by changing 5^ and 5^. We may 
assume for example that the dark pool has a positive inventory at time u G [t, T]. It can 
make buy orders more attractive than sell orders, so to rebalance its position on the stock. 
In particular, by setting 5^ close to the mid-price and 5^ far away from the mid-price (close 
to the best bid price), it encourages buy orders while sell orders are discouraged. At each 
time u G [t, T], we thus consider three options for the inventory management: i) the dark 
poofs order flow may be controlled by optimally choosing the per-share commission size, 
ii) a limit order—of which execution is uncertain—is posted in the lit pool, or hi) a costly 
market order is submitted to the lit pool. Our goal is to obtain the critical levels of the 
inventory for which market orders or limit orders are optimal. 

We fix a filtered probability space satisfying the usual conditions 

and augmented by all P-null sets. We model the LOB bid-ask half spread by a continuous- 
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time Markov chain with discrete state space K := {0,5,25...n5}. The chain 

is generated by {Q} = (r^j) such that P[/c(u + du) = j\k[u) = i] = Hjdu and P[/c(u + du) = 
i\k(u) = i] = 1 + Tjjdu, with > 0 for all j / i and The LOB mid-price is 

defined by 

rh rh 

S(ti)=s + J ju(u,S(u))du + J (j{u,S{u))dW{u), (2.1) 

where is a standard -Brownian motion. At ti G [t, T], the best bid and 

ask prices are given by S^(ti) = S(ti) — /c(ti) and S^(ti) = S(ti) + k{ti), respectively. 

We let a and b be R+ U {0}-valued i.i.d. random variables with finite second moment. 
These variables model the size of incoming sell and buy orders in the dark pool. We define 
the inventory process by 


rt^ 


X(ti) = x+ /(u,X(u),a)diV“(u)- 


/ (u,X(u), b) diV^(u), 


( 2 . 2 ) 


where {iV“(u)} and {iV^(u)} are independent Poisson processes with intensities = A“(5“) 
and Ag = A^(5^), respectively. We allow for short-selling and thus, at any time u e [t, T], 
we have sign[X(u)] = {—1,0,1}. Furthermore, we model the cash process {y(u)}t<u< 7 ’ by 


rh 


nti) = y + 


fi (u,Xiu),Y(u),S(ji), 5^(u), b)diV^(u) 
/i (u, A:(u), y (u), S(u), 5“(u), a) diV'Cu). 


(2.3) 


Equations ( |2.2p and ( |2.3P are strongly connected. For example, an incoming buy or¬ 
der of size b reduces the inventory by f (u,X{u),b) and increases the cash amount by 
fi(u,X{u),Y{u),S{u),5^,b). The analogous holds for incoming sell orders. For technical 
reasons, we consider a bounded domain for both, the inventory and the cash processes, by 
letting (x,y) G ^ c R^.We let ^ := [X ,X] x [T, T], where X, X, 7, and Y are real-valued 
constants. Such an assumption is supported by the following financial interpretation: the 
market maker, i.e. the dark pool, is subject to regulations constraints (e.g. internal risk- 
management) which make it hard to hold or short-sell an amount of shares bigger than a 
fixed authorised quantity. Also, the market maker exits the market whenever the cash pro¬ 
cess has reached either the lower bound (bankruptcy) or an upper bound (gains target). 
The dark pool can resort to the lit pool to liquidate (respectively refill) part of the inventory; 
we assume that it can not post speculative orders. This means that at time u G [t,T] a buy 
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order can be posted only if X(u) < 0 while a sell order can be posted if X(u) > 0. 

Limit orders strategy 

The dark pool (market maker) can post a limit order by specifying a quantity rj and a limit 
price S =b (k + k), where k is the optimal distance from the best price, at which it wants 
to buy or sell. We only consider immediate-or-cancel orders and we model their execution 
percentage by a [0,1]-valued random variable In particular, if a limit order is posted at 
time Tj, for j = 1,2..., then it impacts the inventory and the cash processes as follows: 

Y(Tj)=Y(Tj-)+xiVj,XiTj),Z,S^TjlKTjlKj). 

For a limit buy order, the function F is non-negative valued and the function x is non¬ 
positive valued. In fact, a limit buy order, if executed, increases the inventory and reduces 
the cash amount. The contrary holds for limit sell orders. We will state these assumptions 
rigorously in the appendix. We let be the set of stopping times not greater than T, and Jf 
be the set of all admissible control actions. A limit-order strategy is a collection of stopping 
times and actions L = (tj, rjj, Kj) g x ^ x [0, k]. 

Market orders strategy 

Alternatively, the market maker can submit a market order, which (i) is more expensive and 
(ii) benefits of sure execution as it is matched with existing limit orders. A market order of 
size posted at a time pi impacts the inventory and cash processes as follows 


X(pO=X(p;-)+A(?i,^(pO). ^2 5) 

Y{pi)=Y (pi-)+c(^i,X{pi), S(pi), fc(p;)). 

For a market buy order, the function A is non-negative and the function c is non-positive. 
The contrary holds for market sell orders. A market-order strategy is a collection of stopping 
times and actions M = (pj, x 3^. 

3 Optimisation problem and viscosity solution 

We consider the problem of maximising expected terminal cash subject to total liquidation 
(the problem is the same for a pre-specified non-zero terminal inventory) of the remaining 
inventory via a market order. In defining the objective function, we are led by Guilbaud & 
Pham (2013). We let the stopping time t* be the first time the state variables exit the set 


5 



6, such that 


T* = inf {u > 1 1 (X(u), 7(u)) ^0\KT. 


We define the value function by 


y (^t,x,y,s;k) = sup E 

D,L,M 


g(u,X(u))du 


t<Pi<T 


E 

t<Tj<T 


+ t/(X(T*),7(T*),S(T*),/c(T*)) 


(3.1) 


where D = (5*^, 5^), the function U is the time-T* liquidation revenues and g is a running 
penalty for the risk of holding the inventory In particular, (i) if g is negative-valued, the 
dark pool is risk-averse, (ii) if g is equal to zero, then it is risk-neutral, and (iii) if g is 
positive, then the dark pool is risk-prone. In the summations we include the penalties for 
trading in the lit pool, where > e/ > 1. Throughout the paper we write, whenever 
possible, the vector of state variables x := [x,y,s] G ^ x =: y and we let be the 
set of all stopping times less than t*. Equation ( |3.1| ) satisfies the DPI^ see Fleming & Soner 
(2006). That is, for all t g 


y [t,x;k) = sup E 

D,L,M 


r g(u,x(u))du- X X ^ 1 +^ 

Jt t<pi<T t<Tj<T 


(3.2) 

This is an optimal double-obstacles impulse control problem. We define the non-local oper¬ 
ators if and for limit and market orders respectively, by 


ify 


and 


it,x;k)= sup V(t,x + r(i 7 ,x,z),j + /(i 7 ,x,z,s,/c,K-),s;fc) 4 ''^(z)dz-e;, 

rieJ^,Ke[0,K] Jq 


^V(t,x;k) = supV (t,x + A(^,x),y + c(^,x,s,k),s;k) -e^. 


We define the operator {t,x,k,p,q,M,4)^ = sup^, {t,x,k,p,q,M, (t>,D) by 


j 4 (t,x,k,p,a,M,(t>) 




where 


6 






^(t,x,p,q,M) :=p + ju(t,5)q + |c7^(t,5)M, 

^ 00 

:= sup Ag (0 (t,x+/(t,x,a),j -/i(t,x,5“,a),s;fc) - ^ (t,x; fc)) £a(a)da, 
5 “e[ 0 ,J:] Jo 

r 00 

^l(.t,x,4>):= sup Ag (t,x-/(t,x,b),J+/l(t,x,5^b),s;^:) - ^ (t,x;fc)) £b(b)db, 
5 '’e[ 0 ,fc] Jo 

^4>(t,x;k)= ^ rfcfc/(</)(t,x;fc')-</)(t,x;/c)). 

k'^k 

The value function V{t,x;k) satisfies the QVI system 


min {-g(t,x) - ^ (t,x, /c, d,V, d,V, d,,V, V) XV - JiV){t,x-,k)XV - ^V) (t,x; k)} = 0, 

(3.3) 

on [o,r) X 5^ X K. Equation ( [3.3| ) can be interpreted as follows: if V — > 0 and 

V — 5£V >0, then the value function can not be improved by an impulse and thus no orders 
are submitted to the lit pool. As soon as V — < 0 or V — <0, the value function is 

settoV—= Oory—ify = 0 and an impulse takes place. We thus consider intervention 
times (tj and pj) and impulses (rjj and by which the dark pool can control the evolution 
of the state variables X{u) and Y{u). For this purpose, we define the continuation region 
(CR), the limit orders impulse region (L7) and the market orders impulse region (Ml) by 


CR := {iu,x,k) G [0, T) X ^ X K : y > and y > J^V}, 

LI := {(u, jc, k) G [0, T) X ^ X K : .^y = y and > JIV }, (3.4) 

MI := {{u, x,k)e[0,T)xy xK: MV = y and MV > L£V \. 

For each k G K, let us consider the upper and lower semi-continuous envelopes of the 
function y(*; k) defined by 


y*(t,jc; k) = lim sup y(t^,jc^;k), 


V^it,x;k) = 


lim inf V(t\x^;k). 


Definition 3.1. A system of functions V : [0, r)x^xK^Risa viscosity subsolution, (resp. 
supersolution), of^if 


min git, x)- jV (t,x, k, 8^4), d^cj), dss<p,<p) ,iV*-J(V*) (t, x; fc) , 

iV*-Sev*){i,x;k)^ <0, 


1 





^resp. min|-g(t,x)-^ 5^0,5^0,4^,0) ,[y*-JlV^') {i,x;k) , 

(K-^K=)(t,x;fc)}>oj, 

where 0 G T) x y xK) is such that y*(t,x,/c) — 0(t,x;/c) (resp. V^(t,x,k) — 

(l)(t,x;k)) attains its maximum (resp. minimum) at (i,x,k) (Z [t,T) x y xK. 

Definition 3.2. A system of functions V : [0, T) x y xK^R is a viscosity subsolution (resp. 
supersolution) of^if 

min |- g(t,x)- j^(,i,x, k,p,q,M, V*), (V*-^V*Xi,x;k), fc)| < 0, 

(resp.min |-g(t,x)-^(t,x,fc,p,q,M,yj,(y*-^yj(t,x;fc),(y*-i!f\4)(t,x;fc)| > O), 

where 

(i) the function (p € ^^^>^’^([0,7) x 5^ x K) is such that V*(t,x;k) — p(t,x;k) (resp. 
V^it,x;k) — (l)(t,x;k)) attains its maximum (resp. minimum) at (i,x,k) c [0,T) x 
5^ X K, 


(ii) (p,q,M) € ^^’'^V*(t,x;k) (resp. V^(t,x;k)), where the set is the closure 

of the second-order parabolic superjet ^2-+^ defined by 

^^’^V(i,x;k) := | (p,q,M)eR^ : V (t,x;fc) < V(t,x;fc)+ p(t - t) 

+ q(s-s)+ |m(s-s)^|, 

(Hi) the set S?^~V — —V) is the second-order parabolic subjet. 


We use Definition 3.1 to prove the existence of a viscosity solution and Definition 3.2 


to 


prove the strong comparison result, which implies uniqueness, for discontinuous viscosity 
solutions. The results and the standing assumptions are found in the appendix. We note 
further that, since the control space is bounded, the terminal condition is regular and equal 
to V[T,x;k) = U[x,k). For further details, see e.g. Pham (2009). 
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4 Numerical results 


In this section we explore some explicit examples of the impulse-control problem presented 
in this work and we analyse in details their qualitative features. We start by presenting 
a toy model, for which there exists a semi-explicit solution. Then we move forward by 
progressively removing the simplifying assumptions listed next. For the time being, we 
suppose that 

(a) the mid-price process follows an arithmetic Brownian motion, 

(b) dark-pool orders are unit-sized, 

(c) lit-pool orders are unit-sized, 

(d) the dark-pool commissions 5^ and 5^ are fixed a priori, 

(e) limit orders can not be partially executed (i.e. ^ G {0,1}), 

(f) limit orders can only be posted at the best bid-ask prices (i.e. k = 0), 

(g) the dark pool observes the system and can intervene at fixed discrete times, 

(h) g(u,x) = 0 and U{x,y,s) = y -\-sx — kx^, 

(i) the bid-ask half spread is constant and set to k. 

In assumption (h), g(u, x) = 0 means that the dark pool is risk-neutral. Moreover, given 
the choice of the function U, it seeks to maximise the terminal cash amount, subject to full 
liquidation of its inventory via a market order. Let us fix t = Iq < ti < t 2 < - • < = T 

to be an equally-spaced partition of the time interval T — t, such that tj+i — = A > 0 for 

i = 0,1,...,N — 1. We define the observed inventory and cash processes by 

= YitO - (S(t^+i) - + (S(t;+i) + 5“)(iV“(t;+i) - iV“(t^)). 


In the case that trades are executed only within the dark pool, the value function of the thus 
uncontrolled problem simplifies to 

V(t, X, y,s) = E [y(T) +X(T)S(r) - kX\T)] = y + xs + c-^iT - t)-k^x + 03(7 - t)f, 


9 



where —/c(A^+A^) and C 2 = (A^ —A^). Next, we introduce the possibility 

of submitting unit-sized market and limit orders in the lit pool. At each time tj, the dark 
pool checks whether it is more convenient to (i) execute trades in the dark pool only (ii) 
submit a market order such that 




Y(ti) = y(ti-) - (s(ti) + , 

where = {+ 1 , — 1 } and > 0 , or (hi) submit a limit order such that 




Y(ti) = Y(ti-) - r](ti) (S(ti) - kr](ti)) z, 

where = {+1,-1} and ^ is a {0, l}-valued random variable with P[^ = 1] = £(zi) 
and P[^ = 0] = £(zq). Within this particular example, a different time-ordering of the same 
actions does not influence the value function. We thus provide the following definition. 

Definition 4.1. Let n < N and let Q := ^ be respectively the number 

of market sell (MS), market buy (MB), limit sell (LS) and limit buy (LB) orders submitted in 
[T — (n — 1)A, T], such that < n — 1. We say that two strategies Qi and 

Q 2 are distinguishable iflQi^ Qi and Q 2 ^ Q 2 such that Qi / Q 2 . 

We have n* = ]^^^Q(n+i) distinguishable strategies for the time-interval [T—(n—1)A, T]. 

The situation can be summarised as follows: we have four objects and n — 1 slots. We need 
to count how many combinations we may have, including repetitions, provided that the 
order does not count and some slots may stay empty. Such number n* is readily obtained: 

(n-l) (n-l-i) (n-l-i-j) (n-l-i-j-k) - 3 

z z z 

i=0 j=0 k=0 h=0 i=0 

By the DPI^ we can solve the above problem by means of backward recursion. In particular, 
let Q be the optimal strategy for the interval [T — (n — 1)A, T]. Then the value function is 
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specified by 


y (T - (n - 1)A,X, j,5, Ql,Q^,Q%Q^) = y + xs + c,{T - (n - 1)A) - (/c + e„) (Q^ + Q^) 

+ (kp - ei) (Q“ + Qf) -k{x + (Qf - Q‘l)p + 

-Q^ + C2(T-(n-l)A))", 

whereP[z = 1] = p andP[z = 0] = 1—p. The value function V(T—n A, x,j,s,Q^,Q[|j,QpQ|’) 
l/(r-nA,Q^,Q^,Q“,Q|’) is given by 


LB :l/(Q“„,Q^,Q?,Qf + l) 
DP 

iS :V(Q“„,Q^,Q“ + l,Qf) 
MS:y(Q^ + l,Q^,Q“,Qf) 


■X >, 


if X < min -i —^ — 1 — x, 

\ 2k 

e;-2/c-e^ p _ €i 1-p _ 

if -:-- - X < X < -1-X, 


€i-2k-e^ p 
2fc(l-p) 2 


2fc(l-p) 2 

€n 


2kp 


■c I - e; 1-p _ 

if max t-;— 1 — X, —;--X 

' 2k 2kp 2 

ei 1-p 


< X < min < — + 1 — X, , 

' 2k 2kp 




.r 1-p - -ei + 2/c + e^ p _ 

11 -1-X < X < -tTT^-:-^-X, 


2kp 


.r ^ I + P 


2fc(l-p) 2 

+ ?-^K 


where x = (q|’ — Q“)p + + C 2 (T — nA). Limit buy orders and limit sell orders are 

never the optimal choice if e; < p(e^ + 3/c+pA:). Take for example limit buy orders. At each 
time tj, they are the optimal choice if and only if the inventory x satisfies 


€i-2k-e^ p _ e, 1-p 

-77-1-X < X <-;-- 

2ic(l-p) 2 2kp 2 


(4.1) 


For such X to exist, one must have 


€i-2k- 
2k(l -p) 


-e; + 

2kp 


0 < 6/ < p(e^ + 3/c + pk). 


The same condition also applies for the existence of limit sell orders. Since the above condi¬ 
tion neither depends on time, nor on the strategy followed thus far, we may regard it as an 
intertemporal and strategy-independent condition. Following an analogous procedure, we 
find that the above conditions ensure the existence of limit sell orders, too. In Figure we 
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Value Function 


plot the value function obtained by the recursive relation given above. In Figurewe show 
that the condition e/ < p(e^ + 3k + p/c) ensures that limit orders are the optimal choice 
when the inventory satisfies Equation In particular, this means that limit orders are 

optimal if the penalty for trading in the lit pool by means of limit orders does not exceed the 
quantity p(e^ + 3/c + p/c) as in Figure[^(a). On the contrary in Figure[^(b) we note that an 
optimal combination of market orders and dark pool activity outperforms limit orders and 
thus the latter are never the optimal choice. 


25 -period QVI Optimal Boundary 



Figure 1: n(n + l)(n -F 2)(n + 3)/24 = 20,475 distinguishable combinations of market, limit and dark pool 
orders. We set 5^ = 5^ = 0.25, = 0.5, k = 0.5, = 5, = 3, p = 0.92. 


25 -period QVI Optimal Boundary 


25 -period QVI Optimal Boundary 



Figure 2: Value function for dark pool only (dashed line), for limit orders and dark pool (dotted line), for 
market orders and dark pool (solid line). Limit orders perform better than market orders if < p(e^-F3/c-Fpk). 
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We now remove assumptions (a), (b), (f), (g), and (h) and replace them with 

(a’) the mid-price process follows a geometric Brownian motion, 

(b’) the dark-pool orders a and b are of random size with supports A and B respectively, 

(f) the limit price at which the dark pool posts in the lit pool can be optimally chosen, i.e. 
K^O, 

(g’) the dark pool can observe and control its cash and inventory levels on a continuous¬ 
time basis, 

(h’) g(u,x) = —(pgix) = —(px^ = and U{x,y,s) = y+ x(s —/c|x|/x). That is, the dark pool 
is subject to a quadratic running penalty for holding a non-zero inventory and liquidates 
all its inventory at the final date T. This manifests the dark pool’s risk-adversity since 
a large inventory reduces the value function of the pool. 


We assume that P[^ = 0] = and P[^ = 1] = f= 1 — This reflects the fact 

that the filling-probability of a limit order depends on how far from the mid-price such an 
order is posted. The associated QVI is 


minj(^g(x)-5fy --a^s^dssV 


- J ^V(,t,x- a,y + a(s + 5“),5) - V(t,x,y,s)j ia(a)da 

- J [v(t,x + b, j - bis - 5^),s) - Vit,x,y,s)jlhib)db; 


Vit,x,y,s)- sup [V {t,x + ^,y - ^is + ^k),s) - e^] ; 
?=±i 


Vit,x,y,s)- sup 

ri=±l,KG[0, 


V [V {t,x + r]Zi,y-Zi (s-rjik + K)) ,s) - ei] =0. 


Given the form of the terminal condition, we consider the ansatz y(t,x,y,s) = y + xs + 
/z(t, x) with terminal condition h(T, x) = —/c|x| — e^. We refer to Cartea & Jaimungal (2012) 
and Cartea et al. (2014) for more details about this ansatz. We then get: 
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min |0g(x) — ~ j |^a5“+ /i(t,x — a) —£a(a)da 

[b5^+h(t,x + b) - /i(t,x)] tb(b)db; 




Kt,x)- sup [-\^\k-€^+hit,x + ^)]; 
?=±i 


Kt,x)- sup V [|i 7 |(fc + >c)zi-e,+h(t,x + i 7 Zi)]£''(z;)| =0. 

if]=±l,Ke[0,K] i=Q J 


50-period QVI Optimal Boundary 



■Market Sell Orders 

■ Limit Sell Orders k=0 

■ Limit Sell Orders k=1 

■ Limit Sell Orders k=2 
HDark Pool 

□ Limit Buy Orders k=2 

□ Limit Buy Orders k=1 

■ Limit Buy Orders k=0 


Time 


Figure 3: Optimal boundary for the doube-obstacle problem with k: = 0,1,2. 


As one would expect, as the inventory increases, the dark pool will first resort to limit orders 
and ultimately to market orders. We note that when posting a limit order, the smaller the 
inventory, the farther from the mid-price the market maker posts. A pre-specified non-zero 
inventory would produce a shift in the optimal boundaries by an equal amount. Moreover, 
the symmetry in Figure is due to the particular choice made for the frequency of orders 
submitted by the clients to the dark pool. That is, by choosing / A^, we would loose 
such a symmetry and see a shift upwards (resp. downwards) of the optimal boundaries 
when A^ < A^ (resp. (A^ > A^). To solve the QVI, finite difference methods were applied. 
We now replace assumption (d) with 

(d’) the dark pool optimally chooses the commissions 5^ e {5^, 5^} and 5^ e {5^, 5^}. 
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The QVI now reads: 


min I </)g(x) — 5fh(t,x) — sup A“(5“) | + h(t,x — a) — h(t,x)l£a(a)da 

i 5-=5i Ja*- -■ 

sup A^(5^) I + /i(t,x + b) —/i(t,x)l£f,(^’)db; 

5>’=5l Jb 


5>>=5^ 

b(t,x)- sup [-|?|fc-e„ + b(t,x + <^)] ; 

?=±i 


b(t,x)- sup V [|T]|(fc + K)Zi-e,+b(t,x + i7Zi)] £'"(Zi)| =0. 

r)=±l,Ke[0,k] J 

In Figurewe plot the numerical solution of (4.2). 


50-period QVI Optimal Boundary 



Market Buy Orders 
Limit Buy Orders k=0 
Limit Buy Orders k=1 
Limit Buy Orders k=2 

Dark Pool 5^=5^ 5^=5'^ 
+ + 

Dark Pool 5^=5^ 5'^=5'^ 
+ - 

Dark Pool 5^=5^ 5'^=5^ 

Dark Pool 5^=5^ 5'^=5'^ 
- + 

Dark Pool 5^=5^ 5'^=5'^ 
+ + 

Limit Sell Orders k=2 
Limit Sell Orders k=1 
Limit Sell Orders k=0 
Market Sell Orders 


(4.2) 


20 25 30 

Time 


Figure 4: Optimal boundary for the doube-obstacle problem with k: = 0 ,1,2 and 5°_ = 5^ = 0.2 and = 
5l = 0.4. 

In Figure]^ we see that when the inventory is relatively small, it is optimal to set 5^ = 5^ 
and 5^ = 5^ so to encourage the arrival of dark pool orders. Then, for example, when 
the inventory increases, the incentive is lowered to 5^ = 5^ and 5^ = 5^. If the inventory 
increases further, 5^ = 5^ and 5^ = 5^ is preferred since the dark pool must increase the 
commissions to avoid reaching the threshold at which it will need to place orders in the lit 
pool. The critical inventory level at which the dark pool begins placing orders in the lit pool 
falls as the terminal liquidation date is approached. Finally, we observe that the situation 
is symmetric when the inventory is negative, i.e. the commissions increase as the inventory 
decreases and the lit pool thresholds shrink near the terminal date. In the last simulation, 
we remove assumption (i) and replace it with 

(f) a stochastic bid-ask spread is introduced. 
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For the ansatz V(t, x, y,5; fc) = y + X5 + x), the system of QVIs can be reduced to 


min •! </)g(x) — 5fhj.(t,x) — sup A“(5“) + hj.(t,x — a) — hj.(t,x) fa(a)da 

I 5“=5^ Ja*- ^ 

- sup A^(5^) I rb5^ + hfc(t,x + b)-hfc(t,x)l£b(b)db-Vrfcj rhj(t,x)-hfc(t,x)l; 

5'’=5i Jb'- ^ L J 


hfc(t,x)-sup [-|^|fc-en, + hfc(t,x + ^)] ; 

hfc(t,x)- sup ^ [|i7|(fc + x)Zi-e;+hfc(t,x + i7Zi)]f''(z;)| = 0 . 

?7,Ke[0,/c] J 

(4.3) 


50-period QVI Optimal Boundary 



Time 


50-period QVI Optimal Boundary 



Time 


(a) 1st regime with k = 0.8 


(b) 2nd regime with k = 1.3 


50-period QVI Optimal Boundary 



□ 


I Market Sell Orders 
iLimit Sell Orders k=0 
iLimit Sell Orders k=1 
iLimit Sell Orders k= 2 
-Dark Pool 5^=5^ 5^=5*^ 
I + + 

I Dark Pool 5^=5^ 5*^=5^ 

I Dark Pool 5^=5^ 5'^=5^ 

Dark Pool 5^=5^ 5'^=5'^ 
— + 

„Dark Pool 5^=5^ 5*^=5'^ 

ILimit Buy Orders k=2 
ILimit Buy Orders k=1 
ILimit Buy Orders k=0 
I Market Buy Orders 


(c) 3rd regime with fc = 1.8 


Figure 5: Optimal boundary for the doube-obstacle problem with k: = 0 , 1,2 and 5“ = 5* = 0 .2, 5“ = 5* = 
0.4 and = 0.8,1.3,1.8. 
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In Figure we note that as the spread increases, the lit pool thresholds decrease. In fact, 
the market maker tends to resort to the lit pool earlier in order to avoid holding a large 
inventory at time T, which will indeed be very expensive to dispose of. 

5 Concluding remarks 

In the present work we study an optimal market-making problem faced by a dark pool. 
The dark pool earns the commission fee from executing the trade orders within the pool 
placed by its clients. The market agents, who trade through the dark pools, benefit from 
anonymity and advantageous prices. Throughout the activity the pool faces an inventory 
risk, which can be reduced (i) by controlling the width of the dark pool spread, and (ii) by 
resorting to the lit pool via both, market and limit orders. The dark pool activity is preferred 
since it protects from information leakage. Such a feature is modelled via a fixed penalty 
incurred by the dark pool whenever it submits an order to the lit pool. As confirmed by the 
numerical results, the dark pool will refrain from placing orders in a lit pool as long as the 
size of the inventory is relatively small. Whenever a certain level is exceeded, the dark pool 
resorts to the lit pool by means of limit orders. A limit order is more remunerative though its 
execution is uncertain. The dark pool can choose the limit price; we find that the more the 
inventory grows, the closer to the mid-price the dark pool will post. This is reasonable since 
the filling-intensity of limit orders depends on how far from the mid-price they are posted. 
If the inventory becomes critically large, market orders will be preferred instead, which are 
costly but benefit of sure execution. When the end of the market-making activity approaches 
(which, e.g., might be thought of as the end of the trading day) the market-orders region in 
the lit pool widens while the dark-pool and limit-orders regions in the lit pool diminish. In 
fact, the market maker will incur in a higher penalty by liquidating a large inventory at the 
terminal date. These conclusions are obtained by formulating a double obstacle standard 
stochastic and impulse control problem, and by showing that the associated value function is 
the unique viscosity solution of the corresponding dynamic programming equation. Finally, 
we provide four numerical examples with increasing complexity. The numerical analysis is 
similar to the one carried out by Guilbaud & Pham (2013) and our results concerning the 
lit pool are in line with the ones obtained by them. 
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6 Appendix 

Here state all the standing assumptions which are introduced in the modelling setup pre¬ 
sented in this paper. 

(i) the map jU : [0, T] x ^ R and cr : [0, T] x R_^ ^ R_^ satisfy Lipschitz continuity 

|iu - jU (t2,S2) f + \o- - a (t2,S2) |^ < C (|ti - t2p + |Si -S 2 P) 

for all ti, t 2 € [ 0 , r] and 5^,52 € K+; 

(ii) / : [0, r] X [XjX] xE_,_u{0}^E+U{0} satisfies Lipschitz continuity and the linear 
growth conditions 

\f -f {t2,X2,v)\^ < C (|ti - t2p + |Xi -X 2 P) 

|/ < C (l + |xiP) 

for all tj, t 2 € [ 0 , r], Xi, X 2 e [X,X], and v € M; 

(hi) /i : [0, T] X 5^ X [0, fc] X E+ u {0} ^ M+ U {0} satisfies Lipschitz continuity and the 
linear growth conditions 

|/l (ti,Xi,d,v) -/i (t 2 ,X 2 ,d,v)|^ < C (|ti - t 2 l^ + ||xi -X 2 II 2 ) 

|/i ih,Xi,d,v)\^ < C (l + ||xi|Q 
for all ti, t 2 € [ 0 , r], d e [ 0 , fc] and x^, X 2 e 5^; 

(iv) r : ^ X [X,X] X [0,1] ^ M and / : ^ X [X,X] x [0,1] x M+ x [0 ,k] x K ^ R are 
Lipschitz continuous functions satisfying, for M > 0, the following properties: 

(a) 0 < r(ri,x,z) < M and /(') 7 ,x,z,s,fc, k) < 0 if x < 0, for all 17 e J{, z e [0,1], 
(s, fc) e R_,_ X [ 0 , fc] and fc e K, 

(b) —M < r{r],x,z) < 0 and /(') 7 ,x,z,s,fc, k) > 0 if x > 0 , for all rj e J{, z € [ 0 , 1 ], 
(s, k) € R_,_ X [0, k] and fc € K; 

(c) |xp — E[|x + r(i 7 ,x,z)p] > 1, for all 17 € J{, z € [0,1] and x € [X,X]; 

(v) K ■. SC X [X,X'] —> R and c : SC x [X,X] x R_,_ x K —> R are Lipschitz continuous 
functions satisfying, for M > 0, the following properties: 


18 



(a) 0 < A(i^, x)<M and x,s,k) <0 if x < 0, for all e s e ]R_,_ and fc e K, 

(b) —M < A(i^, x) < 0 and x,s,k)>0 if x > 0, for all e 5 e R+ and fc e K; 

(c) |xp — \x + A(i^,x)p > 1 , for all e ^ and x e [X,A]; 

(vi) g : [0, r] X [X,A] ^ M satisfies Lipschitz continuity 

|g (ti,Xi) - g (t2,X2) 1 ^ < C (|ti - t2p + |xi - X2p) 

for all ti, t 2 € [0, r] and Xi, X 2 € [X ,X]; 

(vii) U : y xK^M. satisfies Lipschitz continuity and the linear growth conditions 

\u{x^,k) - U (x 2 ,/c) I" < C (||xi - X2|0 

|[/(xi,fc)|"<c(l + ||xi||2), 

for all Xi, X 2 ^ y and /c G K. 


6.1 Viscosity Solution |6.1| 

Proposition 6.1. (Existence) The system of functions V(t,x;k) is a viscosity solution of the 
QWQ). 


Proof We use definition 3.1 and we show that the system of functions V(t,x; k) is a viscosity 
solution by proving that it is both a supersolution and a subsolution. First we note that we 
have V(T,x;k) = U(x;k) on {7} x 5^ x K, thus we need to prove the viscosity property 
only on [t,T) x y x K. Results in, e.g., Ly Vath et al. (2007) ensure that 
and i!f\4 < (ifV)^. By definition of the value function, we have V > MV and V >T£V for 
all (u,x) G [t, T) X It follows that V* > {^JIV\ > and V* > > T£V^. That is, 

it suffices to show that 


-g(t,x)-^ (l,x,k,dt(t),ds(t),dss(t>,4>) >0. 

Let (\4 — 0)(t,x;fc) = 0, where [i,x,k) = argmin(\4 — 0)(t,x;/c). By definition of \4, 
there exists a sequence it^,x^) (t,x) such that V^{t^,x^;k) V^(t,x;k) as m ^ oo. 

We define the stopping time 

inf |u > |Xf^j^^(u) ^ B^(t^,x^)j-, (6.1) 


19 





where is the open ball of radius rj centred in We choose a strictly 

positive sequence ^ 0 and let the stopping time 6^ := 0^ A (t^ + A 0* A where 
0 * is the first time the regime switches from its initial value k and where is the first time 
an impulse takes place. By the first part of the DPP and the definition of the function (p, we 
have for any admissible control strategy 




61 




LJl 


An application of Ito’s formula to (p between and 0^ yields 






LJ t 


k, dfCt), d,4), d,,4>, V*, D)] du 


+ . 


We can divide by then let m ^ oo and apply the the mean value theorem. Finally the 
result follows form the arbitrariness of the control variable. First note that if V* < MV"" 
or y* < if y*, the subsolution property is immediately satisfied. We assume therefore that 
y* > .^y* and y* > if y*; we then need to show that 


-g(t,x)-^ (l,x,k,dt4>,ds(l),dss4>,<p) ^0- 


( 6 . 2 ) 


By continuity of the mapping in ( |6.2p , we assume on the contrary that there exist a 5 > 0 
and an e > 0 such that —g[i,x)— j 4 (i,x,k,d^(p,ds(p,dss(p,(p^ > 5, for all G 

B^(^i,x). We take the sequences ^ 0 and ^ valued in B^(^i,x) and we 

define the stopping times 0^ by ( [6.1[ ) with rj < e and 0^ := 0^ A (t^ + h^) A 0* A t^. By 
Ito’s formula and the second part of the DPP, there exists an admissible control strategy D* 
for which 


5K 


r 


Tm 


< E 


[g ^ (u,Xf^ ..^(u), k,dt<P,ds(l),dss<p,(l>,D*)] du 


where = (y* — k). Dividing by —we find that 
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Since 1E[0^ — ^ 1 as m ^ oo, we get 5/2 < 0, which contradicts 5 > 0. □ 


6.2 Strong Comparison Results |6.2| 


Proposition 6.2. (Strong Comparison Principle) Let V and U be a supersolution and a sub¬ 


solution respectively of the QVI (3.3). If U{T,-) <V{T,-\ then U < V on [0,T] x y x 


Proof Let v = v* and u = u* be a supersolution and a subsolution respectively. We first prove 
that there exists a (^-strict supersolution, where 0 < (^ < e;. We refer to, e.g., Seydel (2009) 
for technical details. We consider the function v^(t,x;k) = v(t,x;k) + + |x|^^), 

where /? > 0 and p > 1 are to be determined later. Then we have: 


v^{t,x;k) - ^v^(t,x;/c) 


> v(t,x;fc) + “^^(l + |xpP) — Mv{t,x-,IC) 

- sup (l + |x + hfl, x)P^’) - €^] 




x\^P - sup (|x + A(?,x)pP) 




where the last inequality follows form the supersolution property of the function v. Further¬ 
more, by assumption (v) for the function A, we have 


^eA^-0 


|xpP — sup Hx + A(^,x)pP) 


+ ^171 > 


since |a| > |b| => \af > |b|^Vp > 1. Analogously, we have: 
v^(t,x; fc) — ifv^(t,x; fc) 


> Ce^(^ 


-0 


\x\^P - sup 

k:€[0,?c] Jo 


(|x + r(r],x,z)|2^’)£W(z)dz 


>C. 


Finally we take into consideration the PIDE part. We let be the test function for . Then 
(j) := (j)^ — + |x|^^) is the test function for v. We therefore have: 
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> -g(t,x)-j 2 / (t,x,k,dt(t),ds4>,dss(t>,v) + PCe^^'^~'^\l + \x\'^P) 

r 00 

- sup Ag ^e^^'^~^H\x+f(t,x,a)\'^P-\x\^P)ia(a)da 
5 “e[ 0 ,j:] Jo 

^ 00 

- sup aH (|x-/(t,x, b)pp - |xpp)£i,(b)db, 

5 '’e[ 0 ,J:] Jo 

which implies that 

- g(t,x) - (t,x, fc, dt<p^, dss<p^, v^) 

r 00 

>+ |xpP) — sup Ag Hx+/(t,x,a)pP — |xp^] £a(a)da 

5 “e[ 0 ,J:] Jo 

^ 00 

- sup aH Ce^‘^^"''^(|x-/(t,x, b)p^’-|xpP)£i,(b)db>C, 

5 '’e[ 0 ,J:] Jo 

for [3 sufficiently large. Now we set 


Vm=(l -]v + -vf, u„=[l + -]u-yf. 

\ m m \ m m 


Using Definition 3.1 one can prove that 

min |- g(t,x)- j 4 {t,x, k, 0'") ,(v^-^v^) (t,x; fc), 

y 

(v^-^V;„)(t,x;fc)| > —. 

where 6^ := is the test function for and d) is the test function for v and 

min |- g(t,x)- (t,x, k, 5, (/?'", 5,, (/;'") ,{u^-J(u^){t,x-,k), 

y 

(Um-ifUm)(t,x;fc)} <-, 

J m 

where := is the test function for Ujj^ and Lp is the test function for u. We 

further note that u and v are polynomially bounded (see e.g. Crisafi & Macrina (2014), 
Proposition 6.2, for details). Thus, we have for each /c G K 

lim (u^-v^)(t,x;fc)= lim f 1 +—\u-v)(t,x;k) -+ \x\^P) =-oo, 
x^±oo ^ x^iboo \ mj m 
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where we set p larger than the bounding polynomial of u and v. Thus the supremum is 
attained in a bounded set. Since is upper semicontinuous, it attains a maximum 

over a compact set. Next we show that, for all m large, we have 

M :=max(u^(t,x;k)-v^it,x;k)) < 0, 

t,x,k 

where {i,x,k) = argmax(Un,(t,x;fc) —We define the auxiliary function by 
(ti,t 2 ,Xi,X 2 ;/c) :=u„ (ti,xi;fc) - (t2,X2;k) - T (^|t^ _ t 2 p + ||xi -X 2 II 2 ) , 

For each /c G K, is upper semicontinuous and therefore it admits a maximum 
at ,t 2 ^,x \’^Let be defined by = max^^j^M^’^, attained at the point 

(tp t 2 ,x^^,X 2 ,k^) [i, i,x,x,k) as e ^ 0. Furthermore we have that > M and ^ 

M as e ^ 0. Let us assume on the contrary that > 0. We now go through the various 
cases. Let 

(u^-^u^)(tl,xl;k'^) <0. 

By the supersolution property of and by subtracting the two inequalities, we have 

C 

(u^-^u^)(tl,xl;k^) - - ^v^)(tl,xl;k^) + — <0 

We can now develop a contradiction argument since 
M = lmun,(t^,x^;fc^) - Vn,(t^,x^;fc^) 

c c c 

< lim.^u^rtf,JCT;/c^) — ^v^(t^,xi;k^) -< limM^-= M-. 

e^o V i i y K z z j ^ ^ 

The second case arises when (u^ — /c^) < 0. We follow the same procedure to 

show that 


M = limum{tl,xl;k^) - v^(tl,xl;k^) 


c c c 

< limifu„(tf,xf;fc^) -^v„(d,x^;F)-< limM^-= M-. 

e^o V 1 1 y V z z y ^ ^ ^ 

Next we consider the PIDE part and we use Definition 3.2 Thanks to Crandall & Ishifs 
Lemma, there exist vectors such that 
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Recalling Definition |3.2| we subtract the argument of the subsolution from the argument of 
the supersolution and obtain 

Since we assumed > 0, we choose a g > 0 such that 

0 < qM^ = Q V^(t2,x|;/c^)- ^ (|tj - t2p + ||x^ -X2II2) j 

< g(tl, xl)+J^ ^ 7 (tl - tD’ 7 (^1 - -^1 ’ “m j 

- g(tl,xl) - ^ ^tl,xl,k'^,^(tl - -sl),N^,v^^ 

We can now analyse every component in detail. First we note that, due to Assumption (vi), 

g(tl,xl) - g(tl,xl) < \g(tl,xl)-g(tl,xl)\<c{\tl-tl\ + \xl-xl\) -^0 
as e ^ 0. Furthermore, by Crandall & Ishii’s Lemma and Assumption (i), we have 

- lc72(t^,S^)iV|j < £ + |s^ . 

For the integral part, we analyse . The case involving 

(tpX^,u^) — 3^ [t 2 ,X 2 ,VTn) be treated analogously Given that sup(A) —sup(B) < 
sup(A — B), we have: 


f(tl,xl,u^) - (t|,x^,v„) 

I* 00 

(u^{tl,xl+f{tl,xl,a),yl-fj(tl,xl,5‘^,a),sl;k^)-u^(tl,xl;k^) 

Jo 

- v„(t|,x| +/(t|,x|,a),y| -/i(t|,x^,5“,a),s|;F) +v^(t^,x^;fc^) j£a(a)da- 


< sup Xg 

5“e[0,F] 
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The argument of the integral above can be rewritten as 


t2,Xj +f{tl,xl,a),X2 +f{t2,X2,a),yl -/i(tpXp5“,a), -/i(t2,X2,5“,a), 

s^,5|; fc^) - (tt, t2^x^,x^; fc^) - T Qx,^ - xlf + |y[ - jII' 

- 1^1 + / (t^ a) - x| - / (t^, x|, a) 1^- \yl - (t^,x^, 5“ a) - y| + /i (t^, x^, 5“ a) 

^ ^ (1^1 - ^2r + r - ki + / (ti’a) - ^2 -/ (t2>^1’«) r 

- |jr -/i(t^^^5“,a) - j|+/i(t^,x^,5“,a)|^), 

where the inequality is justified by the fact that the function attains its maximum at 
(tp t|,XpX|;fc^). That is, 


< 


(t^,x^,u,„) - (t^,X^,V,„) 

sup A“ f l-^(\xl-xl\^ + \yl-y^ 2 f-\xl+f(tl,xl,a)-xl-f(tl,x^ 2 ’^)\^ 

“e[ 0 ,fc"] Jo i- ^ 

- \yl - fi (t^ a) - y| +/i (t|,x^, 5“, a) jfaCa) da. 


Since the right-hand-side of the above equality tends to zero as e ^ 0, we have that 
lim^^o “m) “ (^ 2 ’ •’‘^ 2 ’ '^m) ^ 0- Finally, we have 




v^{tl,xl;k^)'j <0 


since the maximum is attained at k^. Thus, by letting e ^ 0, we get qM < 0, which is 
a contradiction since Q > 0. Therefore M < 0. Furthermore, since we have proved that 
u* < v^, the value function is continuous as it is both upper and lower semicontinuous. □ 
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